Fulde-Ferrell superfluidity in ultracold Fermi gases with Rashba spin-orbit coupling 
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We theoretically investigate the inhomogeneous Fulde-Ferrell (FF) superfluidity in a three dimen- 
sional atomic Fermi gas with Rashba spin-orbit coupling near a broad Feshbach resonance. We show 
that within mean-field theory the FF superfluid state is always more favorable than the standard 
Bardeen-Cooper-Schrieffer (BCS) superfluid state when an in-plane Zeeman fleld is applied. We 
present a qualitative flnite-temperature phase diagram near resonance and argue that the predicted 
FF superfluid is observable with experimentally attainable temperatures (i.e., T ~ 0.2Tf, where Tf 
is the characteristic Fermi degenerate temperature). 
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I. INTRODUCTION 

The Fulde-Ferrell (FF) superfluid is a fascinating state 
proposed to understand the fermionic superfluidity with 
unequal populations in the two spin states [1]. Un- 
like the standard Bardeen-Cooper-Schrieffer (BCS) su- 
perfluid, where fermions of opposite spin and momentum 
form Cooper pairs and condense into a microscopic state 
at rest, the FF superfluid is characterized by Cooper pairs 
carrying a single-valued center-of-mass momentum and 
thus by an inhomogeneous condensate state. More com- 
plicated inhomogeneous condensate states are also pos- 
sible with the inclusion of more center-of-mass momenta 
for their spatial structure in real space, as suggested by 
Larkin and Ovchinnikov (LO) [2]. These forms of inho- 
mogeneous superfluidity have now been collectively re- 
ferred to as FFLO superfluids. They are anticipated to 
have manifestations in a number of important physical 
settings, ranging from solid-state superconductors to the 
nuclear matter at the heart of neutron stars tSnSi] . How- 
ever, despite tremendous theoretical and experimental 
efforts over the past fifty years, conclusive experimen- 
tal evidence of their existence remains elusive |3l-[l7|. In 
this work, we show that the long-sought FF superfluid 
might be observable in a three dimensional (3D) ultra- 
cold atomic Fermi gas with Rashba spin-orbit coupling 
and in-plane Zeeman field. 

The idea that inhomogeneous superfluidity is enhanced 
by Rashba spin-orbit coupling was first suggested by 
Barzykin and Gor'kov in the study of surface super- 
fluidity in materials such as WOaiNa [T^, [ly] and was 
later generalized to a 3D Fermi system by Agterberg and 
Kaur [20]. It was recently revisited by Zheng and co- 
workers in the context of ultracold atomic Fermi gases 
|2ll , [22I , which have the unique experimental advantage 
of unprecedented controllability in interactions, spin- 
populations and purity [23]. In contrast to solid-state 
superconductors, ultracold atomic Fermi gases are gen- 
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Figure 1: (color online) Finite temperature phase diagram 
of a 3D Rashba spin-orbit coupled atomic Fermi gas at a 
broad Feshbach resonance {l/kpcis = 0) and at the spin-orbit 
coupling strength Xkr/Ep = 1. The critical temperature, 
in units of the mean-field critical temperature without spin- 
orbit coupling Tco — 0.496Tf, decreases monotonically with 
increasing the in-plane Zeeman field h. 



erally prepared in the strongly interacting regime [5] , the 
so-called crossover regime from a Bose-Einstein conden- 
sate (BEC) to a BCS superfluid, in order to have an 
experimentally attainable superfluid transition tempera- 
ture. As a result, the theoretical investigation of inhomo- 
geneous superfluidity in such ultracold matter has to rely 
on heavy numerical calculations within mean-field the- 
ory. In the previous work by Zheng and co-workers [21], 
the FF superfluidity has been addressed at zero temper- 
ature along the BEC-BCS crossover. Here, we present a 
qualitative mean-field phase diagram at finite tempera- 
tures and show that the FF superfluid state might be 
within reach at the typical experimental temperature 
T ^ 0.2Tf ^ 0.4Tco, where Tf is the Fermi degenerate 
temperature and Tco — 0.496Tf is the mean- field critical 
temperature of a resonantly interacting Fermi gas with- 
out spin-orbit coupling. Our main result is summarized 



in Fig. [H where the phase transition temperature to a 
normal state is shown as a function of the in-plane Zee- 
man field, for a strongly interacting Rashba spin-orbit 
coupled Fermi gas at a broad Feshbach resonance. We 
find that the FF superfluid is always the true ground 
state at a finite Zeeman field. 

It is important to note that the experimentally real- 
ized spin-orbit coupling is not of the pure Rashba type 
[2J,[25|. Instead, it is an equal weight combination of the 
Rashba and Dresselhaus spin-orbit couplings (for a de- 
tailed discussion, see for example, Ref. [26]). The possi- 
bility of observing FF superfluidity in current experimen- 
tal settings has been discussed by Shenoy [27], the present 
authors [28], and also by Wu and co-workers but in two 
dimensions [29]. We anticipate that the Rashba spin- 
orbit coupling might be experimentally realized soon [32] . 
Furthermore, it is also feasible to create a 3D isotropic 
spin-orbit coupling [33]. The FF superfluidity with 3D 
isotropic spin-orbit coupling has been investigated most 
recently by Dong and co-workers [30] and by Zhou and 
colleagues |31l] . 

The reminder of the paper is organized as follows. In 
Sec. II we introduce the model Hamiltonian for a Rashba 
spin-orbit coupled Fermi gas with an in-plane Zeeman 
field and describe the mean-field framework. We present 
the explicit expression of the mean-field thermodynamic 
potential, with the FF pairing momentum q and pair- 
ing order parameter A as the variational parameters. In 
Sec. Ill, we discuss in detail competing ground states 
near a broad Feshbach resonance and show that the FF 
superfluid is always the true ground state at finite Zee- 
man fields in the superfluid phase. We explore system- 
atically the properties of this exotic state of matter at 
finite temperatures. In Sec. IV we present our conclu- 
sions. For the convenience of numerical calculations, we 
use a non-standard form of Rashba spin-orbit coupling. 
In the Appendix A, we show that it is fully equivalent 
to the Rashba spin-orbit coupling commonly used in the 
literature. 



II. MODEL HAMILTONIAN 

We consider a 3D spin-1/2 Fermi gas of ^Li or ^^K 
atoms near broad Feshbach resonances with Rashba- type 
spin-orbit coupling \{(Jxkx -^ (^z^z) and an in-plane Zee- 
man field along the z-direction /la^, a configuration to be 
experimentally realized in the future [32]. Here cix and 
(Jz are the Pauli matrices, kx = —idx and kz = —idz are 
the momentum operators. In the Appendix A, we discuss 
in more detail about the expression of Rashba spin-orbit 
coupling. The model Hamiltonian of the system may be 
described by 



where the single-particle Hamiltonian takes the form 



H= /6/x[Ho(x)+Hint(x)], 
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and the pairing interaction Hamiltonian is given by 



Hh 



C/o4(x)^|(x)^^(x)^t(> 



(2) 



(3) 



describing the contact interaction between the two spin 
states with interaction strength Uq. In the above Hamil- 
tonian, ipl (x) and i/ja- (x) are respectively the creation 
and annihilation field operators for atoms in the spin- 
state cr, and ^k = —h^V^/{2m) — /i is the single-particle 
kinetic energy with atomic mass m in reference to the 
chemical potential ji. We have denoted the strength of 
Rashba spin-orbit coupling and of in-plane Zeeman field 
by A and /i, respectively. The use of the contact inter- 
atomic interaction generally leads to an ultraviolet diver- 
gence at large momentum and high energy. To remove 
such a divergence, it is useful to express the interaction 
strength Uq in terms of the 5-wave scattering length a^ , 



m 



Uq Airh^as V 
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(4) 



where V is the volume of the system. The scattering 
length as can be tuned precisely to arbitrary value, by 
sweeping an external magnetic field across the Feshbach 
resonance l23l. 



A. Mean-field Bogoliubov-de Gennes theory 

A solid-state Fermi system with Rashba spin-orbit cou- 
pling and in-plane Zeeman field provides a promising 
platform to observe the long-sought FFLO superfluid, 
as suggested by Barzykin and Gor'kov in their pioneer- 
ing work [18]. In the context of ultracold atomic Fermi 
gases, this idea was renewed, as motivated by the inter- 
esting finding that the two-body bound state of the model 
Hamiltonian acquires a finite center-of-mass momentum 
along the z-axis ^]. This strongly indicates the exis- 
tence of a FF pairing state with a single- valued center- 
of-mass momentum at the many-body level [23|. Indeed, 
in a recent zero-temperature calculation by Zheng and 
co-workers [21], the parameter space of the FF super- 
fluid has been found to be significantly enhanced by the 
Rashba spin-orbit coupling. Here we explore the whole 
mean- field phase diagram at finite temperatures. 

Let us assume an order parameter with a single- valued 
center-of-mass momentum along the z-axis: 



A(x) = -C/o(^;(x)^tW) = Ae^' 



(5) 



(1) 



The direction of the FF pairing momentum is chosen fol- 
lowing the center-of-mass momentum of the two-particle 
ground state [34]. It is also consistent with the previous 
mean-field studies for Rashba spin-orbit coupled Fermi 



systems [ISl, [20|, [III- Within mean- field theory, we ap- 
proximate the interaction Hamiltonian by, 



-Hi: 



A(x)^|(x)V'|(x)+H.c. 



i^- <^' 



For a Fermi superfluid, it is convenient to use the fohow- 
ing Nambu spinor representation for the field operators. 



$(x)= Ut(x),^;(x),^|(x),^j(x)' 



(7) 



where the first two and last two field operators in ^(x) 
could be interpreted as the annihilation operators for par- 



J 



tides and holes, respectively. The total Hamiltonian can 
then be written in a compact form. 



n = ^l dx$t (x)?^BrfG*(x) - ^^ + E ?'k, (8) 



where the factor of 1/2 in the first term arises from the 
double use of particle and hole operators in the Nambu 
spinor <^(x). Accordingly, a zero-point energy X^^^k ap- 
pears in the last term, which is formally divergent. The 
Bogoliubov Hamiltonian I-LbcLG takes the form. 



^BdG 



& -^Xkz^h \kx -A (x) 

\kx i\^-\kz-h A (x) 

A* (x) -4 -^Xk^-h \k^ 

-A*(x) A^^ -i\.-\k^h\ 



(9) 



In free space, where the momentum is a good quantum number, it is straightforward to diagonalize the Bogoliubov 
Hamiltonian 



'^BdG^kTyCx) = £;kr7^kr7(x), 

by using the plane-wave quasiparticle wave-function 



V 



(10) 



(11) 



and quasiparticle energy E\^r)' The Bogoliubov Hamiltonian now becomes a 4 by 4 matrix and the four eigenvalues 
and eigenstates have been specified using the index 77 (77 = 1, 2, 3, 4). 

The mean-field thermodynamic potential 1] at a temperature T is then given by 
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^ (^k+q/2 + ^k-q/2) - ^ El^V 
k k77 



keT 
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k77 
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where the zero-point energy Xl^^ E\^r] in the first term 
(i.e., the square bracket) is again due to the double use 
of particle and hole operators, and the zero-point en- 
ergy Ek^k has been rewritten as Ek(^k+q/2 + ^k-q/2)/2 
to cancel the divergence in Ekrv^kry- The second term 
in the above thermodynamic potential accounts for the 
thermal excitations of Bogoliubov quasiparticles, which 
do not interact with each other in the mean-field approx- 
imation. It should be noted that, the summation over 
the quasiparticle energy in Ekry ^^^st be restricted to 
E\ir] ^ 0, because of an inherent particle- hole symme- 
try in the Nambu spinor representation. For instance, 
it is straightforward to check that for any particle state 
[^kt,^ki,^kt,^ki]^ with energy Ey, > 0, there is a one- 
to-one corresponding hole state ['^1^^, ^^k^' '^-kf '^-k^]^ 
with energy —E\^. These two states correspond to the 



same physical solution. 

It is easy to show that, in the absence of spin-orbit cou- 
pling (A = 0), we may explicitly write down the expres- 
sion for the quasiparticle energy E\^r)- Eq. (p!2|) then re- 
covers the thermodynamic potential of a spin-imbalanced 
3D Fermi gas 0. 



B. Numerical solutions 

For a given set of parameters (i.e., the temperature T, 
the interaction parameter l/kpCLs etc.), we solve the or- 
der parameter A and the FF pairing momentum q by us- 
ing the self-consistent stationary conditions: dfl/dA = 
and dVt/dq = 0, together with the number equation 
N = —dQ/d/i for the chemical potential /i. Different 
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Figure 2: (color online) The free energy of three competing 
states, including (i) a normal Fermi gas with A = 0; (ii) a 
fully paired BCS superfluid with A 7^ and q = 0; and (iii) a 
finite- momentum paired FF superfluid with A 7^ and ^ 7^ 0, 
as a function of the in-plane Zeeman field at resonance and 
at T = 0.05Tf. The Rashba spin-orbit coupling strength is 
Xkp /Ef = 1. The inset shows the Zeeman field dependence 
of the BCS (red solid line) and FF order parameters (black 
solid circles). 



saddle-point solutions of these coupled equations give 
competing ground states. At finite temperatures, the 
true ground state is the one that has the lowest free en- 
ergy F = ^ + (iN. 

In numerical calculations, we take the Fermi wave- 
vector kp and the Fermi energy Ep as the units for 
wave-vector and energy, respectively, in order to make 
the equations dimensionless. We focus on the unitary 
limit with a divergent scattering length l/ikpCis) = 0. 
Throughout the paper, we shall use a Rashba spin-orbit 
coupling strength Xkr/Ep = 1. We will consider the 
superfluid transition temperature as a function of the in- 
plane Zeeman field, as well as the critical Zeeman field 
across the Feshbach resonance at a given temperature. 



III. RESULTS AND DISCUSSIONS 

For any set of parameters, in general there are three 
competing states: the normal gas (A = 0), BCS super- 
fluid (A 7^ and g = 0), and FF superfluid (A 7^ and 
q ^ 0). These competing states all satisfy the stability 
condition d^Vt/d/S? > and therefore are stable against 
phase separation in real space. In Fig. [21 we show the free 
energy of these states in the unitary limit at T = O.OSTp. 
It is readily seen that the FF superfluid is always more 
favorable in energy than the standard BCS superfluid at 
a finite in-plane Zeeman field, when the Rashba spin- 
orbit coupling is present. This observation holds for any 
interaction parameters and temperatures. 

At small in-plane Zeeman fields (i.e., h < 0.5Ef), the 
FF superfluid and BCS superfluid are very similar, with 



essentially the same pairing gap, as shown in the inset of 
Fig. [21 The difference in the free energy between these 
two superfluids is indeed negligible. However, strictly 
speaking, we always find dVt/dq < for the BCS su- 
perfluid, indicating that a FF superfluid with a finite 
center-of-mass momentum q is preferable, although the 
magnitude of q might be very small. With increasing the 
Zeeman field, the BCS superfluid ceases to exist. The FF 
superfluid will also disappear, but at a bit larger critical 
Zeeman field. The FF pairing gap vanishes much more 
smoothly than the BCS pairing gap, as can be seen from 
the inset of Fig. [2l The phase transition from the FF 
superfluid to the normal state is continuous. 

The significantly enlarged parameter space for a FF su- 
perfluid can be qualitatively understood from the change 
of the two Fermi surfaces due to spin-orbit coupling, as 
discussed by Barzykin and Gor'kov [18| and by Agterberg 
and Kaur [20]. Here we explain very briefly the physical 
picture following their ideas. With our choice of Rashba 
spin-orbit coupling, the single-particle energy spectrum 
takes the form, 



E\i± 



2m 



±^\^kl^{\k,^hf, 



(13) 



where "±" accounts for two helicity branches. At rela- 
tively large spin-orbit coupling (i.e., Xkp ^ h)^ we may 
approximate 



E^ 



h^kl 



k± 



where 



2m 



^(^k^±rj'±xVk[Tki, (14) 



2m.h 



h^,/W+kl, 



(15) 



It is easy to see that the centers of the two Fermi surfaces 
are shifted by q/2 along the z-axis in opposite directions. 
In this situation, we may shift the two Fermi surfaces 
by an amount ±{q/2)ez^ respectively [20|, to maximize 
the overlap between the two surfaces. Here e^ is the 
unit vector along the z-axis. The fermionic pairing then 
occurs between the single-particle states of k + {q/2)ez 
and — k + ((7/2)6^, giving rise to a FF order parameter 
that has a spatial dependence A(x) = Ae**^^. This pair- 
ing mechanism holds for arbitrary small in-plane Zeeman 
field and explains why the FF superfluid is always more 
preferable than the standard BCS superfluid. The direc- 
tion of the FF momentum is uniquely determined by the 
form of spin-orbit coupling and its magnitude is roughly 
proportional to the in-plane Zeeman field /i, i.e., see Eq. 

By collecting the critical in-plane Zeeman field for the 
FF superfluity state at different finite temperatures, we 
obtain a finite-temperature phase diagram in the strongly 
interacting unitary limit, as shown in Fig. [H For a 
2D weakly interacting Fermi gas with Rashba spin-orbit 
coupling, such a phase diagram has been determined by 
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Figure 3: (color online) The free energy difference between 
the FF superfluid state and the possible BCS superfluid state 
at resonance and at three different temperatures, T = 0.05Tf 
(black solid line), O.IOTf (red dashed line), and 0.20Tf (blue 
dot-dashed line). The Rashba spin-orbit coupling strength is 
Xkp / Ef — 1- The BCS state ceases to exist above a critical 
Zeeman field he — 0.72Ef, as indicated by the thin vertical 
line. At h > he, we replace Fbcs by Fn in the calculation of 
the free energy difference. 
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Figure 4: (color online) The pairing momentum of the FF 
superfluid at resonance and at three different temperatures, 
T = O.OSTf (black solid line), O.IOTf (red dashed line), and 
0.20Tf (blue dot-dashed line). The Rashba spin-orbit cou- 
pling strength is Xkp /Ep = 1. 



A. Temperature dependence of the FF superfluid 



Barzykin and Gor'kov [18]. At small Zeeman fields /i ~ 0, 
we find that 



[T, ih = 0)- TJ /T, {h = 0)cx h\ (16) 

in agreement with the analytic expression in Ref. [18]. At 
the low-temperature regime, the transition temperature 
Tc{h) is roughly given by 



T,(/i)a[/ie(T = 0)-/i]", 



(17) 



where a > 2 and hc(T = 0) is the critical Zeeman field 
at zero temperature. Our result of a > 2 seems to be 
consistent with the analytic prediction of a = 3 for a 2D 
Rashba spin-orbit coupled Fermi gas |l8l |. 

It is important to note that our phase diagram in Fig. 
[1] seems to be qualitatively different from the one ob- 
tained in the previous study by Zheng and collaborators, 
which predicts a BCS superfluid at small in-plane Zee- 
man field [21]. In addition, the critical Zeeman field at 
zero temperature hcJT = 0) is shown in their Fig. 1(d) to 
be about 0.55Ef [2F|, much smaller than what we have 
obtained numerically, i.e., hc{T = 0) ~ Ejr. The quali- 
tatively different phase diagram is simply due to the dif- 
ferent interpretation of the FF state with small pairing 
momentum. For example, in the calculations by Zheng 
and collaborators, the FF state with pairing momentum 

M: 10~^kF has been regarded as the BCS superfluid 
. On the other hand, the different critical Zeeman 
field is probably due to the different accuracy of numer- 
ical calculations 13511. 



We now explore in greater detail the finite-temperature 
properties of the FF superfluid. Fig. [3] reports the differ- 
ence in free energy between the FF superfluid state and 
the BCS superfluid state (the normal state) at three typi- 
cal temperatures ^]. The energy difference per particle 
is sizable at low temperatures, suggesting that the FF 
superfluid is very robust with respect to other compet- 
ing ground states. The difference becomes considerably 
smaller with increasing temperature. Nevertheless, it is 
still visible (i.e., at about 0.002Ef) at the typical exper- 
imental temperature T = 0.2Tf. Thus, the thermody- 
namic stability of the FF superfluid could be guaranteed. 

Fig. m shows the FF pairing momentum as a function 
of the in-plane Zeeman field at different temperatures. At 
low Zeeman fields, the pairing momentum is essentially 
independent on the temperature. At high fields close to 
the critical Zeeman field, however, the pairing momen- 
tum decreases quickly with increasing temperature. We 
note that, near zero temperature, the maximum pairing 
momentum is comparable to the Fermi wave- vector kp- 

Fig. [5] presents the population imbalance between the 
two spin states, calculated by using V5n = —dVt/dh and 
normalized by the density np = A:|./(37r^). The BCS 
superfluid state has low capacity to accommodate the 
population imbalance. Thus, it ceases to exist at a criti- 
cal Zeeman field. The deformation of the Fermi surfaces 
in the FF superfluid state is able to allow more popu- 
lation imbalances. At low temperatures (upper panel, 
T = O.OSTp), we observe that the population imbalance 
of the BCS superfluid rises quickly near the critical Zee- 
man field, compared with that of the FF superfluid. This 
distinct behavior might be used to identify the FF su- 
perfluid. At a relatively large temperature (low panel. 
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Figure 5: (color online) The number difference between the 
two spin states at T = 0.05Tf (upper panel) and T — 0.20Tf 
(lower panel) as a function of the in-plane Zeeman field at 
resonance and at the spin-orbit coupling strength Xkr/Er = 
1. In each panel, the number difference of the three competing 
states is shown by solid circles (FF), red line (BCS), and blue 
dot-dashed line (normal Fermi gas). 



T = 0.2Tf), the difference in the population imbalance 
between the two superfluids is smeared out by the tem- 
perature. 



B. FF superfluid across the BEC-BCS crossover 

We have so far focused on the resonance limit with 
l/{kFCis) = 0. In Fig. [6l we show the critical Zeeman 
field for the FF superfluid state across the BEC-BCS 
crossover at the temperature T = O.IT^. It increases 
monotonically when the Fermi cloud crosses from the 
BCS limit to the BEC limit. In the inset, we present 
the FF pairing momentum at the critical Zeeman field. 
There is a maximum in the pairing momentum at about 
the Feshbach resonance, due to the competition between 
the many-body and two-body effects. The initial increase 
of the FF pairing momentum on the BCS side arises from 
the many-body effect. Towards the BEC limit, however, 
the two-body effect becomes dominant and the FF pair- 
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Figure 6: (color online) Phase diagram along the BEC-BCS 
crossover at T = 0.1 OTf and Xkp /Ef — 1- The main figure 
reports the critical Zeeman field at which the FF superfluid 
phase turns into the normal state. The inset shows the FF 
pairing momentum at the critical Zeeman field. This figure 
might be contrasted with the zero temperature result in Ref. 
[21], see for example, Fig. 1(c) in Ref. |21l |. 



ing momentum decreases gradually to the center-of-mass 
momentum of the two-particle bound state. In the pre- 
vious study by Dong and collaborators, the many-body 
and two-body predictions for the FF pairing momentum 
have been compared with each other for a 3D Fermi gas 
with 3D isotropic spin-orbit coupling [30]. Our result in 
the inset is consistent with their predictions. 



C. Phase diagram at smaU spin-orbit coupling 

Before closing this section, we briefly discuss the phase 
diagram at small spin-orbit coupling. In the absence of 
spin-orbit coupling, the formation of an inhomogeneous 
superfluid is driven by the population imbalance, which 
also leads to the distortion of the Fermi surfaces. In that 
case, there are many equivalent ways to deform the sur- 
faces. As a result, the direction of the single FF pairing 
momentum is not specified. Thus, for any FF super- 
fiuid solution with a pairing momentum +q, we can al- 
ways find another degenerate solution with the pairing 
momentum — q. This indicates that a stripe LO phase 
with an order parameter in the form of cos(q • x) will be 
more favorable, which is simply a superposition of the 
+q and — q plane waves. The investigation of the LO 
phase in a 3D Fermi gas without spin-orbit coupling has 
been carried out by Burkhardt and Rainer many years 
ago [37] and recently by a number of authors [38, 39]. 
In the presence of spin-orbit coupling, the two solutions 
with the +q and — q plane waves are no longer degen- 
erate. Numerically, we find that one solution becomes 
more favorable and the energy difference between the 
two solutions increases rapidly with increasing spin-orbit 



coupling. Therefore, with increasing spin-orbit couphng, 
we anticipate the LO phase wih cease to exist and the 
large spin-orbit coupling will uniquely determine a single- 
valued pairing momentum and lead to a FF superfluid. 

We note that the competition between LO and FF 
phases in a Rashba spin-orbit coupled 3D Fermi^gas was 
recently investigated by Agterberg and Kaur [20]. The 
stripe LO phase was found to cease to exist at large pop- 
ulation imbalance. 



IV. CONCLUSIONS 

In summary, we have investigated the finite- 
temperature phase diagram of a three dimensional atomic 
Fermi gas with Rashba spin-orbit coupling and in-plane 
Zeeman field near a broad Feshbach resonance. We have 
shown that its superfiuid state is always an inhomoge- 
neous Fulde-Ferrell superfiuid. 

Our work complements the previous studies with spin- 
orbit coupled Fermi systems in the solid-state, which are 
described by essentially the same model Hamiltonian in 
the weakly interacting regime [l8l-[2Ql|. Our mean-field 
treatment gives a qualitative picture of the Fulde-Ferrell 
superfiuidity in the strongly interacting resonance limit, 
which is of great interest and of experimental relevance. 

Our result extends the previous zero-temperature in- 
vestigation by Zheng and collaborators [21] to finite 
temperatures. In particular, we have clarified that the 
BCS superfiuid found by these authors at small in-plane 
Zeeman fields is better understood as the Fulde-Ferrell 
superfiuid with small pairing momentum. Our finite- 
temperature calculations indicate that the Fulde-Ferrell 
superfiuidity could be observable at about T ~ 0.2Ti?, a 
temperature that is already reached in current cold-atom 
experiments when the spin-orbit coupling is absent [36]. 

Our investigation is based on the mean-field theory, 
which is known to provide qualitative picture of the BEC- 
BCS crossover. The strong pair fluctuations in a strongly 
interacting Fulde-Ferrel superfiuid might be taken into 
account by using many-body T-matrix theories [4Q| - [42| . 
This is to be addressed in the future work. 
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Appendix A: Rashba spin-orbit Hamiltonian 

In the literature, the Rashba spin-orbit coupling and 
in-plane Zeeman field is commonly written as X{aykx — 
(^xky) — h(Jx [21]. Here we show that our single-particle 
Hamiltonian Eq. ([2]) uses exactly the same Rashba-type 
spin-orbit coupling and in-plane Zeeman field, after some 
rotations in real space or spin space. 

For this purpose, we first perform the rotation of real 
space along the x-axis: x ^ x, y ^ — z, and z ^ y. 
As the result, the single-particle Hamiltonian Eq. (|2j) 
becomes. 



^V^l 



Xkx i\^- Xky - h _ 






. (Al) 



In the second step, let us take a unitary transformation. 



1 



V^t (x) = -^ [<^t (x) - h (x)] , 



(A2) 



V^;(x) = — H)[(/>^(x) + (/>;(x)]. (A3) 

It is straightforward to check that under such a transfor- 
mation, the single-particle Hamiltonian changes to. 



/)|,( 



^k —iXkx — Xky — h (j)^ 

iXkx — Xky — h ^k J L ^i 

which precisely has the form of X{(Tykx — CFxky) — hcTx- 

The interaction Hamiltonian Eq. (j3]) is apparently 
not affected by the spatial rotation. Furthermore, as 
V^l (x) ijj^ (x) = —i(j)^ (x) (j)^ (x), it is also not affected by 
the unitary transformation. For the FF pairing momen- 
tum, we can see that the direction along the z-axis is 
fully equivalent to the direction of ?/-axis after the two 
rotations. Thus, our assumed 2:-axis for the FF pairing 
momentum is consistent with the observation by Zheng 
and collaborators [21] and by Barzykin and Gor'kov fl8| , 
that the FF pairing momentum with X{aykx—(Jxky) — hax 
is along the ?/-axis. 
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